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We calculate via the effective field theory (EFT) approach the next-to-next-to-leading order 
(NNLO) spinl-spin2 conservative potential for a binary. Hereby, we first demonstrate the abil- 
ity of the EFT approach to go at NNLO in post-Newtonian (PN) corrections from spin effects. The 
NNLO spinl-spin2 interaction is evaluated at fourth PN order for a binary of maximally rotating 
compact objects. This sector includes contributions from diagrams, which are not pure spinl-spin2 
diagrams, as they contribute through the leading-order spin accelerations and precessions, that 
should be first taken into account here. The fact that the spin is derivative-coupled adds signif- 
icantly to the complexity of computations. In particular, for the irreducible two-loop diagrams, 
. . which are the most complicated to evaluate in this sector, irreducible two-loop tensor integrals up 

Cn , to order 4 are required. The EFT calculation is carried out in terms of the nonrelativistic gravita- 

tional (NRG) fields. However, not all of the benefits of the NRG fields apply to spin interactions, 
as all possible diagram topologies are realized at each order of G included. Still, the NRG fields 
remain advantageous, and thus there was no use of automated computations in this work. Our final 
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I. INTRODUCTION 

Though predictions of General Relativity (GR) have been confirmed in all observations and experiments to date, 
one of its most essential predictions, namely, that of gravitational waves (GW), has not yet been directly observed. 
Worldwide efforts are undertaken in order to observe such signals either with ground-based [iHil , or future space-based 
[5| detectors. Inspiralling compact binaries, which are promising candidate sources for such signals, can be described 
analytically by the post-Newtonian (PN) approximation of GR [0|. It turns out that even relatively high-order PN 
corrections for the binary inspiral, as high as the fourth order PN (4PN) correction, have a phenomenological impact 
on the theoretical waveform templates, required for successful detection (see [3] and references therein). Moreover, 
astrophysical objects are expected to have significant spins, for which spin effects have a big impact on the event rates 
expected in GW detectors [8|. It is therefore desirable to have PN corrections involving spin effects to the same high 
orders as for the nonspinning case. 

A novel effective field theory (EFT) approach was suggested recently by Goldberger and Rothstein for the treatment 
of the binary inspiral problem [9|, |lO[ . The EFT approach is very advantageous in applying the efficient standard 
tools of quantum field theory to GR, such as Feynman diagrams and dimensional regularization. Subsequently, 
PN corrections of conservative dynamics have been reproduced: IPN and 2PN corrections for binaries [3, [ll|, 2PN 
correction for the n-body problem Xl\ , first promoting the use of automated computations in the EFT a ppr oach, and 
recently even the 3PN correction [7|. The EFT approach was further extended to include spin effects [l^, and PN 
corrections involving spin were tackled as well: the next-to-leading (NLO ) sp inl-spin2 |14| - KL7j (complete results also 
in [1^) and NLO spin-squared interactions [T^ (complete results also in |20l- |22| ) at 3PN order were computed. The 
more complex NLO spin-orbit interaction at 2.5PN order was also computed |23l - [25| (first obtained in [26'J29']). We 
note that following [30, l3l| in more traditional methods, the NLO spin-orbit and spin(a)-spin(b) interactions for the 
n-body problem [S^] and the next-to-next-to-leading order spin-orbit interaction for a binary [33 were also obtained. 
Following the EFT treatment of the radiation sector in the nonspinning case fsij, the spin-orbit, spinl-spin2, and 
spin-squared components of multipole moments were computed to NLO in i35] , where spin-orbit radiative effects at 
NLO were already obtained in |23], and have been pushed further to include tail effects in [Sg]. Recently, the EFT 
formalism has been extended to incorporate radiation reaction [37|, and similar EFT approaches were developed to 
investigate the extreme mass ratio inspiral problem [38-40], which is also relevant for GW detection, and to treat 
weak ultra relativistic scattering |41| . 

A major improvement for the obtainment of higher-order PN corrections via the EFT approach was presented in 
|42| . There, a reduction over the time dimension of the metric a la Kaluza-Klein was made, and was demonstrated 
to improve the IPN order EFT computation. This nonrelativistic parametrization of the metric defines a set of 
new nonrelativistic gravitational (NRG) fields, which were later used to reproduce, e.g. the 2PN and 3PN order 
corrections via EFT as was already noted. The advantages of the NRG fields are numerous. First, there is the 
physical interpretation of the different field components and the clear coupling hierarchy to the mass and spin. 
Second, the derivation of the self-gravitational vertices is simple since there is a full explicit expression for the pure 
gravitational action [43|, as well as a simple expression for its stationary part which contributes at leading orders. 
Further, simple propagators are obtained, there are no mixed 2-point functions, and the derivation of the mass 
couplings is also simple and immediate. Because of the structure of vertices and mass couplings to the worldline, not 
all possible diagram topologies are realized at each order, and the number of diagrams, and, in particular, the more 
complicated ones, is reduced. Those are pushed to higher orders, for examplCjWith respect to the standard Lorentz 



covariant parametrization, where a one- loop diagram is eliminated at IPN 42|, or to the Arnowitt-Deser-Misner 
(ADM) parametrization, where mass couplings are eliminated at 2PN |44| . 

In this paper, we calculate the next-to-next-to-leading order (NNLO) spinl-spin2 conservative potential for a binary 
of compact spinning objects at the 4PN order, via the EFT approach in terms of the NRG fields. With this result, we 
demonstrate for the first time the ability of the EFT approach to go beyond the NLO in PN corrections involving spin. 
The NNLO spinl-spin2 interaction sector includes contributions from 51 diagrams, of which 45 are pure spinl-spin2 
diagrams, while further 6 arise from other sectors, that contribute through the LO spin equations of motion (EOM), 
that should first be taken into account here. Of the pure spinl-spin2 diagrams, there are 39 new diagrams here, 
while 6 others already appeared at the NLO spinl-spin2 sector, though they include new ingredients in the worldline 
couplings. In particular, there are seven two-loop diagrams contributing. We note that, unfortunately, not all of the 
benefits of the NRG fields apply for spin interactions, as all possible diagram topologies are realized at each order of G 
included, which was already illustrated in the NLO spin interactions. Still, the NRG fields remain advantageous, and 
thus there was no use of automated computations in this work. Our final result can be reduced, and a corresponding 
Hamiltonian may be derived. An alternative derivation of the NNLO spinl-spin2 Hamiltonian can be found in [4^ . 

Throughout this paper, we use c = 1, r/^^ = (iiag[l, — 1, — 1, — 1], and the convention for the Riemann tensor is 
R'^^j^p = ^qFJ^o — dpV^^ + F^^r^^ — F^oF^Q,. Greek letters denote indices in the global coordinate frame, while 
lowercase Latin letters from the beginning of the alphabet denote indices in the local Lorentz frame. All indices run 



from to 3, while spatial tensor indices from 1 to 3, are denoted with lowercase Latin letters from the middle of the 
alphabet. The notation Jj^ = J rf;;^ is used for abbreviation (boldface characters denote d-dimensional vectors). The 

scalar triple product appears here with no brackets, i.e. ax b- c= (ax b) ■ c (as there is in fact no ambiguity regarding 
the order in which the product can be performed). 

The paper is organized as follows. In Sec.|lTl we briefly review the EFT approach for the binary inspiral with spinning 
objects, and present the Feynman rules required for the EFT computation with the NRG fields. In Sec. lIIIi we present 
the evaluation of the NNLO spinl-spin2 Lagrangian/Routhian, going over all contributing Feynman diagrams, and 
giving the value of each diagram. In Sec. IIVI we present the NNLO spinl-spin2 Lagrangian/Routhian EFT result, 
and explain how to derive from it the NNLO spinl-spin2 Hamiltonian. In Scc.|Vl we summarize our main conclusions. 
Finally, in Appendix |3 we include the tensor Fourier, one-loop and two- loop reduced integrals required here, whereas 
in Appendix [Bl we give the LO spin EOM that contribute here. 

II. EFT APPROACH FOR BINARY INSPIRAL WITH SPINNING OBJECTS 

In this section, we present the ingredients required in order to perform the EFT calculation of the NNLO spinl-spin2 
interaction in terms of NRG fields, namely, the Feynman rules and the effective action from which they are derived. 
Here, we review briefly and build on Sees. II and III of [25| and references therein, following similar notations and 
conventions as those that were used there. 

First, we parametrize the metric in a nonrelativistic form according to the Kaluza-Klein ansatz 

dT^ = g^^dx^'dx" = e^-^idt - A, dx'f ~ e-^'^'j^jdx'dx^ , (1) 

defining the set of nonrelativistic gravitational (NRG) fields (p, Ai, jij = 6ij + aij. Then, in terms of the NRG fields 
the metric reads 

-e'^'t'Aj \^( 1 + 20 + 202 -Aj-2Aj(t)-2AjP 



3m«^ (^ -e^M, -e-2^7,j + e'^A.Aj )-\ -A, - 2A,(j> - 2Ajcf>^ -% + 20% - a,, - 2^^S,j + 20a,, + A,A, J ' 

(2) 
where we have written the approximation for the metric in the weak-field limit up to the orders in the fields that are 
required for this work. 

The action describing the dynamics of the binary system is given by 

D ^= jg + Opp, [o) 

where Sg is the pure gravitational action, and Spp is the worldline point-particle action for each of the two particles 
in the binary. 

We consider first the purely gravitational action. It is the usual Einstein-Hilbert action plus a gauge-fixing term, 
which we take as the fully harmonic gauge, i.e. 

Sg = Seh + Sgf = -Y^/^'^V5^+ 5^/^'^V55p.^^^^ (4) 

where F'^ = T'i^g'"^ . The full explicit expression for the Einstein-Hilbert action and for the fully harmonic gauge 
fixing in terms of NRG fields was given in [4J|, where it was obtained using Cartan's method of 2-forms. Thus, there 
is no need to expand for each required ingredient specifically, and the propagators and self-gravitational vertices can 
be obtained readily from the action. 

The NRG scalar, vector, and 2-tensor field propagators in the harmonic gauge are then given by 

= ( 0(a;i) 0(a;2) ) = 4^^ ^ ""' j^'^ "' ^(ii - h), (5) 

=(A,(xi) v4,(x2)) = -16^G(S.y / ^ ^^ 6iti-t2), (6) 

^^^= = {cT,,{xi)akiix2)) = 327rG P.^^ f "" ^^ '^(^i " ^2), (7) 

where Pij-ki = h i^ik^ji + SuSjk — 2Sijdki)- Here and henceforth, the Feynman rules are presented in position space. 
This makes more sense considering the nature of the binary inspiral problem, in which the external positions of the 
particles are given, rather than the usual external momenta in quantum field theory. 



There are also time-dependent quadratic vertices, which result from the fact that the propagators are actually 
relativistic rather than instantaneous. The Weyl rescaling present in the NRG parametrization eliminates undesired 
mixed quadratic vertices, so that the 2-point functions between the three different fields are 0: {(pAi) — {(p ajk) = 
{AifTjk) = 0. Thus, the Feynman rules for the propagator correction vertices are given by 

^ = 8^ /'^'^ ^^*'^^'' ^^^ 

--X-^-3^/d^x(a.A.f, (9) 

) C = Y2^ /^'^ [^{dt^vf - (dtau)'] , (10) 

where the crosses represent the self-gravitational quadratic vertices, which contain two time derivatives. 

There are also contributions from three-graviton vertices of cubic gravitational self-interaction. The Feynman rules 
for the three-graviton vertices required for the NNLO of the spinl-spin2 interaction are given by 
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^Jd^^^ [d,A, id,A, - d,Ai) + (d^^f) , (11) 

-^ = --^^ I d^^ [^<^^3 {^^Ak^JAk + dkAAAj - 2dkA,djAk + 2^^AJ^kAk) 

(d^Aj {d,Aj - djA,) + i^,A^f)] , (12) 

dx {2(Tijdi(j)dj(f> — ajjdi(f)di(j)) , (13) 



1 



16ttG 




-^ f d^x AAcf,dtcf>, (14) 
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d^x [2cry (diCpdtAj - dt(j>diAj) - ajj {di(pdtAi - dtcpdiAi)] , (15) 

d'^x A.diAjdtAj. (16) 




where the first three vertices are stationary, and can be read off the stationary Kaluza-Klein part of the EH action. 
The three last vertices are time dependent, and each contain a single time derivative. 

In the NNLO spinl-spin2 interaction, contributions from four-graviton vertices of quartic gravitational self- 
interaction also appear. The Feynman rule for the four-graviton vertex required to the order considered here is 
given by 

= ^Jd^^ 0' {d,A, {d,A, - d,A,) + {d,A,f) , (17) 

where this vertex is again stationary. 

Now, we consider the worldline point-particle action. This is given by 

i>pp — ^pp(m) + '-'pp(S)i (18) 

where Spptm) denotes the coupling of the point particles to gravity without the inclusion of spin, and Sppts\ the 
coupling of the spin degrees of freedom of the particles to gravity. Considering the gravitational coupling to the two 
massive compact objects, we take the worldline action of a point particle for each of the objects, so that we have 

2 



Spp{m) = - V'7^i / dXi, (19) 



where finite-size effects are not taken into account here as their contribution enters at higher orders p . We parametrize 
the worldhne using the coordinate time t = a;", i.e. A = i, so that we have for u^ = dx^ /dX: u'^ = 1, u^ = dx'' /dt = w\ 
Thus, the Feynman rules for the one-graviton couphngs to the worldhne mass required for the NNLO of the spinl-spin2 
interaction are given by 



dt (p 



l + -v^ 



m 
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dt Ay (1 + ---) 



dt aijV^v^ (1 
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(20) 



(21) 



(22) 



where the heavy solid lines represent the worldlines and the spherical black blobs represent the particles masses on 
the worldhne. The ellipsis denotes higher orders in v, beyond the order considered here. 

For the two-graviton couplings to the worldhne mass required at this order, we have the following Feynman rule: 



m 



dt02 (1 + ---) 



(23) 



As expected, mass couplings do not play a major role in the spinl-spin2 interaction. 

Next, we consider the gravitational coupling to the two spinning compact objects. Here, we are not concerned with 
finite-size effects of spin, which are quadratic in the individual spins, and hence are not relevant for the spinl-spin2 
interaction (see [19| for spinl-spinl effects). Thus we consider here only the part of the point particle action, which 
is linear in the spin of the particles given by 



5pp(S) 



4E I dKSfcO^abU^ 



(24) 



where 



26 — 



s'^i^,, 



are the Ricci rotation coefficients, and this form for the couplings was introduced in [l] 



following the Routhian in 46]. Since the spin degrees of freedom are naturally formulated in terms of tetrads, here it 
is more convenient to start from the standard Lorentz covariant parametrization t/^i, = rj^i, + h^^, and the derivation 
of the couplings to the worldhne spins is not so immediate as that of the couplings to the masses. Our background 
reference tetrad, expanded in terms of /i^^jy, is given by 



e"- ^ 5°- ^ -/i° h'^hP -\ h°-h''h'' 



(25) 



Using this tetrad, and expanding up to the third order in /i^j^, which is the order required for the NNLO spinl-spin2 
interaction, we obtain the following Lagrangian: 



L 



pp(s) — T^'S"' ha^.bu^ -\- -S°' h^ ( -ha,y. 



/2 I i^uii.a '^aj.i.u ^' 



+-8^"^, h^hi^^ + Kh^^ + -K (/.;;_, - v,1 ^^^ + 



(26) 



Now we should transform from the Lorentz covariant parametrization to the NRG fields using Eq. 1^ in order to 
obtain the couplings to the worldhne spin. 



The Feynman rules for the one-graviton couphngs to the worldhne spin are thus given by 



dt - (S'^d^Aj - S°'d,AjV^ + S^'doA,) 



(27) 



dt {S'WjCJyv' + S°'d,(t) - S°'dQ(t>v') 



(28) 




dt - {S'^d,a,kv'' + S°'doa,y) , 



(29) 



where the (gray) oval blobs represent the spins on the worldlines. Note that here the full expressions for the one- 
graviton spin couplings should be considered. 

For two-graviton couplings to the worldline spin, the Feynman rules required for the NNLO spinl-spin2 interaction 
are: 





dt 



2S'^d^Aj(t) + ]-S'^Ajd,(t) + ]-S'^ Ajdo(t)v' 



+-S°' {Ajdj<l)v' - iAjd.ctiv^ - A.djcjjv^ + djAi<j)v^ - 2d,Aj<j)v^ + A,5o0 + 35o^i 



dt-S'i{djAu-dkAj)a^k + --- 



I dt i^'^" (.AkdjA.v'' + ]^A,d,Auv^ + ^A.duA^v'^ + -^A.d^A^ + is^M, {d,A, - d.Ai) 



(30) 

(31) 

(32) 





dt (25°\ 



dt - [S'^a.k (djcj^v'' + dkc^v^) - S^^'a^.d^cj) +•••], 



(33) 



(34) 



where the ellipsis denotes higher orders in v, beyond the order considered here. The first coupling here was already 
encountered at the NLO spinl-spin2 interaction, but as it is extended here to a higher FN order, it becomes much 
more complicated. Further, new couplings arise here, which are also complex. 

At the NNLO spinl-spin2 interaction, we also have to include three-graviton spin couplings, that are encountered 
for the first time. For three-graviton couplings to the worldline spin, the Feynman rule required here is 




= / dt 4S'*^(/)2a,Aj + -S'^(j)di(l)Aj 



(35) 



where naturally the gravito-magnetic vector is involved in this high-order field coupling LO contribution. 

Two important features of the spin couplings can be noticed very early on and should be stressed. First is the 
fact that the spin, formally being a tensor, is derivative-coupled unlike the scalar mass. This fact translates into 
the requirement of higher-order tensor expressions for all the integrals involved in the calculations (see Appendix |X] 
for more details), which adds significantly to the complexity of computations. Moreover, the derivative-coupling also 



allows for time derivatives in the worldline couplings, which are an additional complication in the spin computations. 
The time derivatives also make the corresponding terms scale at higher PN orders. The second notable feature is the 
fact that the spin couplings contain S"° entries, which represent the redundant unphysical degrees of freedom related 
with the spin tensor. These are taken at this stage as independent degrees of freedom, yet eventually, possibly even 
after the obtainment of the EOM, they are reduced from the final result using some spin supplementary condition 
(SSC). These will also yield contributions of higher PN orders with respect to the S''^ spin tensor components. Both 
features make the PN order of the spin couplings implicit. This makes the power counting, which is essential in the 
EFT approach, more difficult with respect to the nonspinning case. 



III. SPIN1-SPIN2 INTERACTION AT FOURTH POST-NEWTONIAN ORDER 

In this section, we evaluate the relevant two-body effective action by its diagrammatic expansion. For the nPN 
order in spin interactions, we need to consider Feynman diagrams up to the G^'^~^'^ order, where \n~\ is the ceiling 
value of n. Thus in the NNLO spinl-spin2 potential, which is evaluated at 4PN, we have diagram contributions up 
to order G^, coming from all 12 possible topologies appearing at these orders, as displayed in Figs. 1-3 of [llj: one 
topology at 0{G), two at 0{G^), and nine topologies at 0{G^). Unfortunately, not all benefits of the NRG fields are 
present in spin interactions, where in general all topologies are realized at each new order of G included, unlike the 
nonspinning case, where there is a reduction in the number of topologies and diagrams. Hence, using the NRG fields, 
the NLO spin interactions, for example, include the one-loop diagrams, which are omitted from the IPN potential. 
Similarly, the NNLO spinl-spin2 here includes all G^ topologies, such as those with a single cubic vertex or a quartic 



vertex, unlike the NNLO nonspinning case - the 2PN potential computed in |11|. 

For the construction of the Feynman diagrams, we use the Feynman rules from Sec. |TT1 which we PN expand, 
see Sees. IV and V of [2^ for more detail. All in all, we have 51 diagrams contributing to the NNLO spinl-spin2 
interaction, 45 of which are pure spinl-spin2 diagrams. Six other diagrams contribute through the LO EOM from spin 
interactions, as both of the objects are considered spinning: three appeared as spin-orbit diagrams at NLO, and three 
are orbital interaction diagrams, which appeared at 2PN. Of the pure spinl-spin2 diagrams, 39 are new diagrams, 
while 6 others already appeared at the NLO spinl-spin2 sector, although they include new terms in the worldline 
couplings. 

Here we denote r = xi{t) — X2{t), r = \r\, and ft = -. The spin is represented by a 3- vector defined by S^^ = e^^^S^ . 
The labels 1 and 2 are used for the left and right worldlines, respectively. The [1 <->■ 2] notation stands for a similar 
term, whose value is obtained under the interchange of particles labels. Note that under this exchange, n — > — n. 
Finally, a multiplicative factor of J dt is suppressed and omitted from all diagram values. 



A. Order G Feynman diagrams 

For the NNLO spinl-spin2 interaction, we have 13 one-graviton exchange diagrams to evaluate as shown in Figs.[T] 
and [21 Figure [I] contains seven pure spinl-spin2 diagrams, whereas Fig. [5] contains diagrams which appeared in the 
NLO spin-orbit sector or the 2PN orbital interaction sector. Diagrams (a) and (b) in Fig. [l] appeared already in the 
NLO spinl-spin2 evaluation, and they correspond to diagrams (b) and (c), respectively, of Fig. 2 in [la]. The NLO 
evaluation yielded spin-precession terms which were then omitted, but do contribute at this order at the substitution 
of the LO spin-orbit precession, see Appendix [Bl We note that the substitution of lower-order EOM in higher-order 
PN Lagrangians and Hamiltonians is a well-founded procedure, see, e.g. |47H49| . In diagram (b), double-precession 
terms arise, but these will only contribute from the next PN order, so they are dropped here. In principle, new 
diagrams are added in this sector by just inserting further propagator correction vertices. Thus there are three new 
diagrams (as diagrams (d) and (f) here also correspond to (al) and (a2), respectively, of Fig. 2 in 16]). We note that 
unlike the diagrams of nonspinning interactions, which require a tensor Fourier integral of order 2n for n propagator 
correction vertices, here a tensor of order 2?i -I- 2 is required, due to the derivative-coupling of spin, which makes the 
computations heavier. In particular, we note diagram (c), which requires Fourier tensor integrals of orders 5 and 6, 
see, e.g. Appendix VK\ for these Fourier integral tensors up to order 4. The expressions for the orders 5 and 6 are too 
lengthy to be included here, as they have 26 and 76 generic terms, respectively. Again, new precession terms which 
arise at the evaluation at this order are dropped at the use of the LO EOM of spin. We stress that at the NNLO 
spinl-spin2 level, accelerations, and precession terms are inevitable. Finally, we recall that there are several ways to 
evaluate the diagrams including time derivatives, differing by just total time derivatives. 
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FIG. 1. NNLO spinl-spin2 Feynman diagrams of order G: One-graviton exchange. The sohd, dashed, and double hnes represent 
the 4>, Ai, and aij fields, respectively. 



The values of the one-graviton exchange diagrams are then given in the following: 



Fig. 1(a) 
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Fig. l{b) = ^ 51"* ( (52 X a2Y - (52 x n)*a2 ■ n] + S^' { (5i x a^y - (5i x n)'ai • n 
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'?* ((52 X viyd2 • n + (52 X nyd2 ■ vi + {S2 x a2)Vi • n + (52 x ai)V2 • n 
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— 3(5i X V2yvi ■ ftv2 ■ ft + 3(5i X fJ2y{vi ■ nf' + 6(5i x i/i)Vi • ftv2 ■ ft + 6(5i x n) Vi • '!72'S'i • n 
+ 3(5i X ftyv2 ■ ftvi — 3w25i X V2 ■ ftvi ■ ft — 'iv\Si x V2 ■ ftv2 ■ ft + 3n*5i x vi ■ 1/21/2 • ft 

— 3n*5i X V2 ■ ftvi ■ V2 — 15(5i x n)*i72 ■ ^(^1 ■ n) + 15n*5i x V2 ■ ftvi ■ ftv2 -ft) , (37) 



Fig. 1(c) 
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81 ■ 0,282 ■ nai ■ n + 81 ■ n82 • 0,10,2 • n + 81 • nS'2 • 0201 ■ n + 81 ■ 7182 ■ noi ■ 0,2 
35*1 • nS'2 • ndi ■ na,2 ■ n 



G r 

"8^ 



35*1 • 5'2a2 • nv\ + &81 ■ 820,2 ■ wii^i • n — 25*1 • V182 ■ via,2 ■ n — 25"! • V182 ■ 02^1 • n 

— 25*1 • V182 ■ na,2 ■ vi — 28i ■ nS2 ■ via,2 ■ vi — 81 ■ nS'2 • 02^1 — 251 • a2>5'2 • wiwi • n 

— 81 ■ 3,282 ■ nv\ — 98i ■ 82{vi ■ 71)^0,2 ■ n + 681 ■ V182 ■ nvi ■ na,2 ■ n 
+ 65*1 • nS'2 • wi'i^'i • na,2 ■ n + 38i ■ 0,282 ■ n(wi • nf' + 3Si • nS2 • a,2{vi ■ n)^ 
+ 3Si • nS2 • na2 • nvi + 6S1 • nS2 • na2 • i^i^i • n — lb8i ■ 1282 ■ n{vi ■ n)^a2 • n 

— 3Si • ^201 • nv2 — 681 ■ 820,1 ■ W2i/2 • n + 2Si • V282 ■ V2ai ■ n + 28i ■ 0182 ■ W2'i'2 • n 
+ 28i ■ n82 ■ V2ai ■ V2 + 28i ■ V282 ■ noi ■ V2 + 81 ■ 0182 ■ nvl + 28i ■ V2S2 ■ aiV2 ■ n 
+ 81 ■ n82 ■ diV2 + 98i ■ 82{v2 ■ n)^ai • n — 681 ■ 7182 ■ ^2^2 • "^ai • n 

— 681 ■ V282 ■ nv2 ■ ndi ■ n — 3S'i • nS2 • fli (^2 ' nf' — 3Si -0,182 ■ n{v2 ■ nf 

— 3Si • nS2 • ndi ■ nv2 — 6S1 • nS2 • noi ■ "i^ji^z • n + \^8i ■ 7182 ■ n(w2 • n) Oi ■ n 

Q V ^ ^ ^^ ^^ ^^ 

—J 3Si • S2W1W2 + 6S'i • 82{vi ■ V2f - 2S'i • W1S2 • vivl - 45i • V182 ■ W2W1 • V2 

— 4Si • W2S'2 • 1/1 i/i • V2 - 2S'i • V282 ■ V2vl - 98i ■ S2(wi • nfvl - 98i ■ 82{v2 ■ nfvl 

— 36Si • 82V1 ■ 1/2^1 • nv2 ■ n + 6S1 • wiS'2 • wi(i/2 • n) + 12Si • V182 ■ W2W1 • nv2 ■ n 
+ 128i ■ V2S2 ■ 'S'lWi • nw2 • n + 6S1 • W2'S'2 • W2('i''i • nf + 681 ■ V1S2 ■ nvi ■ nw| 

+ 12Si • V1S2 ■ nvi ■ W2W2 • n + 12Si • W2S'2 • nvi ■ W2W1 • n + 6S'i • V282 ■ nw2 • nv\ 
+ 6S1 • nS2 • wi'yi • nw| + 12Si • nS'2 • 'i''i'!7i • W2W2 • n + 125*1 • ^S2 • W2'i''i • 'S2i?i • n 
+ 6S1 • nS2 • W2'i'2 • nv\ + 3Si • nS2 • nv\v2 + 6S1 • nS2 • n(wi • V2f + 45Si • 82{vi ■ nf{v2 ■ nf 

— 30Si • V1S2 ■ nvi ■ n{v2 ■ n) — 305*1 • ^252 • nv2 ■ n{vi ■ n) — 305i • nS2 • wiwi • n{v2 ■ n) 

— 305i • n52 • i72'S2 • ^{^1 ■ nf — 155i • nS2 • n{vi ■ nfv^ — 15Si • n52 • n{v2 ■ nfv\ 

— 6OS1 • n52 • nvi ■ W2W1 • nv2 ■ n + 105Si • nS2 • n{vi ■ nf{v2 ■ nf , (38) 



Fig. 1(d) = ^ 



G 



F^9. 1(e) = g 



8'l'v\S2 X n • 02 - S°*i;^Si x n • ai - 8fv\Sl^n^ + 8°'viS°fn^ 

8^'8°^ (2v\vi - iv\n'v2 ■ n - in^v^vi ■ n) 
— 'i8fv\S2 X V2 ■ nv2 ■ n — 352*W2'5i x ^1 • nvi ■ n , 
81 ■ 820,1 ■ 0,2 — 81 ■ 0,282 ■ oi — 81 X n ■ 0,182 X n ■ 02 



(39) 



AOi cOi 



■ 8'('8^' - S'i"*n*S'2"^n^' + 8^' {{82 x 02)* + n'S2 x n ■ 02) + 8^' ((Si x aif + n'Si x n ■ di 



G 

'2r2 



Si • S2a2 • wit/i • n — Si • S2ai • 'y2'S2 • n + Si • W252 • aiW2 ■ n — 81 ■ 3282 ■ wiwi • n 

— 81 X Vi ■ 7182 X Vi ■ d2 + 81 X V2 ■ di82 X V2 ■ n 

+ 3Si X vi ■ 7182 X n ■ 02 wi • n — 3Si x n ■ di82 x V2 ■ nw2 • n 

+ Sf* (S2"Vi • n + vlS°^n^ + n'S^^vl - 3n'S°^n^vi • n + (S2 x 02)^1 • n - v\S2 x n • 03 

— n*52 X wi • a2 + 3n'S2 x n • 02^1 • n j — S"* ( (S2 x U2) V2 • n + V2S2 x V2 ■ n 

-in' 82 X V2 ■ 7iv2 ■ n] - 8f [sfv2 ■ n + v^Sl^n^ + n' S'l^ v^ - 3n'S"'n^V2 ■ ti 

+ {8i X ai)V2 • n — v\8i x n ■ di — n'Si x V2 ■ di + 3n*5i x n • aiW2 • n ^ 

+5* ( (Si X wi) Vi • n + vlSi X vi ■ n — 3n*5i x vi ■ nvi ■ n ] 
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G_ 

'2^ 



Si ■ 5*2 (wi • ^2)^ - Si ■ V2S2 ■ wii^i ■ V2 + Si X vi ■ V2S2 X ■i/2 • i/i - 35*1 • S2V1 ■ ^/2i/i • nv2 ■ n 

+ 35i • V2S2 ■ wiwi • 'nv2 ■ n — 'SSi x vi ■ V2S2 x V2 ■ nvi ■ n — 3Si x vi ■ nS'2 x V2 • W1W2 • n 

— 35i X vi ■ nS2 X V2 ■ nvi ■ V2 + \hSi x vi ■ 1182 x V2 ■ nvi ■ nv2 ■ n 
+ S'l' {{82 X 1/2)^1 • V2 + ^;2'5'2 X U2 ■ ^^1 - 3(5*2 X U2)Vi • nv2 ■ n - iv\S2 X V2 ■ nv2 ■ n 

— 3w2<S'2 X V2 ■ nvi ■ ft — 3n^S2 x 172 • 1^11/2 • n — 3n^S2 x V2 ■ nvi ■ V2 + 15n*S'2 x V2 ■ nvi ■ nv2 ■ n 
+ S'2* f 01 X viYvi ■ V2 + vlSi X u 1 • 1/2 — 3(51 X ui)Vi • nv2 ■ n — iv^Si x vi ■ nvi ■ n 

— Sv\Si X vi ■ nv2 ■ n ~ 3n^Si x vi ■ ^2^1 • n — 3n^Si x vi ■ nvi ■ V2 
+ 15n*5i X vi ■ nvi ■ nv2 ■ n] 



(40) 



G 



Fig. 1(/) = — 5°^ ((^2 X nya2 ■ vi - 2v{S2 x n • a2 - 4^2 x vi ■ n\ 



G r 



F^9■l{9)-^Yr 



G_ 

'2^ 



82^ ( {Si X nySi ■ V2 — 2v\Si X ft ■ Si — alSi x V2 ■ n 

Si' ({82 X V2yvi ■ V2 + VIS2 X i;2 • Wl 

— 3(52 X n)Vi • W2'U2 • ft — &v\S2 X fj2 ■ fiv2 ■ ft + 3^252 X fji ■ ftv2 ■ ft 
+ 8^' ({Si X ViYvi ■ V2 + vlSi X Ui • W2 

— 3(5i X n) Vi • fl2fJi ■ ft — 6V2S1 X fJi ■ fivi ■ ft + 3w5.5i x V2 ■ ftvi ■ ftj 

Si ■ 0,182 ■ 0,2 — 25i • 0252 ■ 0.1 — 81 ■ ftS2 ■ fioi ■ 0,2 
28i X ft ■ 0,182 X ft ■ CL2 + Si X ft ■ 0,282 X ft ■ 0,1 

Si ■ V182 ■ a2Vi -ft + Si- V1S2 ■ fta2 ■ vi + 81 ■ nS'2 • vid2 ■ vi — 2Si ■ 0,282 ■ wiwi • ft 



(41) 



25l • Wl5^2 • fj2f^l ■ V2~ Si ■ V282 ■ VlfJl ■ V2 + Si X Vi ■ fj2S2 X V2- fJl 



— 25*1 X vi ■ ftS2 X vi ■ 0,2 + 81 X vi ■ 0,282 X vi ■ ft — 35i • nS'2 • na2 • wi^i -ft 
+ QSi X vi ■ ftS2 X ft ■ a2Vi ■ ft — 35*1 x ft ■ 0,282 x vi ■ ftvi ■ ft 

— 81 ■ 0,182 ■ 172^2 ■ ft- Si- fiS2 ■ V2a,i • W2 - 5i • V2S2 ■ ftai ■ V2 + 25i • U2S'2 • a,ifj2 ■ ft 

+ 25*1 X V2 ■ 0,182 X V2 ■ ft — 81 X V2 ■ ftS2 X U2 • fli + 35i • nS'2 • ftdi ■ W2W2 ■ ft 

— 681 X ft ■ 0,182 X V2 ■ ftv2 ■ ft + 35*1 X V2 ■ ft82 X ft ■ diV2 ■ ft 

G_ 

'2^ 

— 38i ■ V1S2 ■ t/2'S'l • ftv2 ■ ft + 65i • U25'2 • t/ll/l • ftv2 ■ ft — 38i • Ui52 • ftf^l ■ t/2'l''2 ■ ft 

— 38i ■ V2S2 ■ ftvi ■ U2i7i • ft — 35i • n5'2 • wi'S'i • i/2W2 • ft — 35*1 • n52 • W2t?i ■ ^2'i''i ■ ft 

— iSi ■ ftS2 ■ ft{vi ■ 1/2)^ — 35*1 X fji ■ V282 X fJi ■ ftfj2 ■ ft — &81 X fji ■ V2S2 X V2 ■ ftfji ■ ft 

— 65i X Vi ■ ftS2 X V2 ■ Vifl2 ■ ft — 35*1 X W2 ■ ^5*2 X V2 ■ fJlfll ■ ft — 681 X di ■ n82 X V2 ■ ddi ■ V2 

+ 35*1 X V2 ■ nS2 X vi ■ nvi ■ V2 + I58i ■ nS2 ■ nvi ■ i/2i7i • nv2 ■ n 

+ 3O81 X vi ■ nS'2 X V2 ■ nvi ■ nv2 ■ n — 155i x V2 ■ nS2 x vi ■ nvi ■ nv2 ■ n 



(42) 



The diagrams appearing in Fig. [2] already appeared in the NLO spin-orbit interaction and the 2PN orbital interac- 
tion. Diagrams (al), (a2), and (b), correspond to diagrams (a3), (bl), and (b2), respectively, in Fig. 2 of j2S]. They 
yield acceleration and precession terms. It can be easily seen, that on substitution of the LO spin-orbit accelerations 
and LO spinl-spin2 precessions, see Appendix iBl NNLO spinl-spin2 interaction terms are obtained. Diagrams (c), 
(dl), and (d2) in Fig. [5] also correspond to diagrams (e), (b), and (c), respectively, in Fig. 4 of [Il|. Again one can 
easily get convinced, that since both objects are considered to be spinning here, on substitution of the LO spinl-spin2 
accelerations, see Appendix [Bl NNLO spinl-spin2 interaction terms are obtained. 
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0--X- 



(y-^^ 



-X- 



-O <> — X — <» <> XX o 



(al) 



(a2) 



(b) 



(c) 



(dl) 



(d2) 



FIG. 2. Feynman diagrams of order G that contribute to the NNLO spinl-spin2 interaction: One-graviton exchange of spin- 
orbit and orbital interactions, which yield acceleration and precession terms. Since both objects are considered to be spinning 
here, these diagrams contribute through the substitution of the LO spin-orbit and spinl-spin2 EOM, see Appendix iBl Diagrams 
(al), (a2), and (b) should be included together with their mirror images. 

The values of the one-graviton exchange diagrams, which contribute through the substitution of EOM, are then 
given by: 



Fig. 2{al) = 2 

r 

Gm2 






Fig. 2(a2) 

Fig. 2{b) 
Fig. 2(c) 



5*1 X ui • a2 — 5*1 X n • a2Vi ■ H ~ Si x ij2 ■ nv2 ■ ft] + Gm2Si y. n ■ 0,2 



—- — ( 5*1 X U2 • fli — 5*1 X n • aiV2 ■ n — Si x vi ■ V2 + Si x vi ■ nv2 ■ n — Si'^vX + 5'?*nV2 • n) 
2r \ I 

— 2G?7iiTO2 (fli • W2W2 ■ n — (12 ■ wi'^i • fi) + 2Gmim2rai ■ 02, 



Fig. 2{dl) — —Gmim2 (ai ■ W1W2 ■ fi ~ 0,2 ■ 1/2^1 • ft) , 

Fig. 2{d2) — —Graim2 (2ai • i'2'i/2 ■ ft + ai ■ ft (^2 — {v2 ' ^) ) ^ 2a2 • fliVi ■ ft — 0,2 ■ ft {v'l — {vi ■ ft) j) 



(43) 

(44) 

(45) 
(46) 
(47) 



- —Gmim2r (ai ■ 0,2 + ai ■ na2 • fi) 
8 



(48) 



As we noted already, there are several ways to evaluate diagrams including time derivatives, as those appearing here, 
e.g. from the NLO spin-orbit interaction. Though these different evaluations lead to physically equivalent potentials, 
one must consistently use the same NLO spin-orbit potential corresponding to a specific evaluation in all stages of 
the calculation, e.g. in the derivation of a corresponding NNLO spinl-spin2 Hamiltonian, where the NLO spin-orbit 
potential should be taken into account. 



B. Order G^ Feynman diagrams 

For the NNLO spinl-spin2 interaction, we have 25 diagrams at order G^ to evaluate: 9 two-graviton exchange 
diagrams, and 16 cubic self-gravitational interaction diagrams as shown in Figs. |3] and SI respectively, (all pure spinl- 
spin2 diagrams). Here, only diagram (a) of Fig. |3] and diagram (al) of Fig. |4] appeared already in the NLO spinl-spin2 
evaluation: they correspond to diagrams (a) and (b), respectively, of Fig. 3 in 16]. All the rest are new diagrams. 

The two-graviton exchange diagrams just factorize into a product of two tensor Fourier integrals, see Appendix VK\ 
We first encounter here two-graviton exchange diagrams, which involve time derivatives, either from the spin couplings 
or from propagator correction vertices. Precession terms that arise here are delegated to higher orders. 

The values of the two-graviton exchange diagrams are given as follows: 



Fig. 3(a) 



+ - 



G2to2 

2r4 
G^m2 



155*1 • S2V2 — 2Si ■ S2V1 ■ ftv2 ■ ft + 2Si ■ ftS2 ■ fJiV2 ■ ft — 395i • nS'2 • ftv2 



Sf 



Fig. 3(6) = 2^ 



?i°'' (3(^2 X viY - 3(52 X V2)' - 2(52 X n)Vi • rt + 10(52 x n)V2 • n - 9n'52 x vi ■ nj 
(5{Si X vi)^ - 5(5^1 X V2)' - 20(5i x n)Vi • rt + 7{Si x n) V2 • rt + Un'Si x V2 ■ ri) , (49) 

f ((5i X V2y - In'Si X V2 ■ «)] , (50) 



Si ■ S2V2 — Si ■ V2S2 ■ V2 — 75i X 772 • ^ 52 X 1/2 • r? + 5' 
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(> 



(g) 

FIG. 3. NNLO spinl-spin2 Feynman diagrams of order G^: Two-graviton exchange. These diagrams should be included 
together with their mirror images. 




Fig. 3(c) 



G^ma r 



2r4 



Si ■ S2V2 + 4:Si ■ S2V1 • i/2 + 5*1 • V1S2 ■ V2- Si- V2S2 ■ V2 + US'! ■ S2{v2 ■ n) 

165i • S2V1 ■ nv2 ■ ft — 135i • V1S2 ■ nv2 ■ n + liSi ■ V2S2 ■ nv2 ■ n — 135i • nS'2 • ■i7iW2 ■ n 
ASi ■ nS'2 • i/2Wi ■ n + ISS*! • nS2 ■ W2W2 ■ n — liSi ■ nS'2 • nvi ■ V2 
63Si ■ nS'2 ■ n{v2 • n)^ + 785i • nS2 ■ nvi ■ nv2 ■ n , 



Fig. 3(d) = - 



G^m2 



2r4 



Si ■ S2V2 ~ 45*1 • S2V1 ■ V2 + 45*1 • V1S2 ■ V2 + 55*1 • V2S2 ■ vi — Si ■ V2S2 ■ V2 — 55i • 52(^2 

165i • S2V1 ■ nv2 ■ n — 135i • V1S2 ■ nv2 ■ n — 205*1 • V2S2 ■ nvi ■ n + 7Si ■ V2S2 ■ nv2 ■ n 
I'SSi ■ nS'2 • 'yiW2 • n — I6S1 ■ nS'2 • 'C2W1 • n + 5Si • nS2 • W2'y2 • n — 13Si • nS2 ■ nvi ■ V2 
15Si • nS2 • n{v2 ■ n)^ + 78Si • nS2 • nwi • nv2 ■ n 



Fig. 3(e) = 2 



G^m2 



4S'i • S2'y| — 2Si • W2S'2 • V2 — 3Si ■ S2{v2 ■ nf + 3S'i • nS2 • ^2^2 • n — 6S1 • nS2 • nw| 



3S'i X V2 ■ n S2 'X V2 ■ n 



Fig. 3(/) - -2 



Fig. 3(5) = -2 



G^m2 



Si ■ S2V2 — Si- §2 {v2 • n)^ + Si • nS2 ■ V2V2 ■ n — Si ■ nS2 ■ nw| 



— 3S'i X V2 ■ n S2 X V2 ■ rt 
/^2 _^ ^^ ^^ ^^ ^^ 

— 5i • 521^2 ^ 2S'i • S'2Wl • W2 - Si • UiS'2 • 1/2 - Si • 'i72S2 • V2 - 3S'i • S2{V2 ■ nY 

+ 5i • S'2Wi • ni72 • n + 3Si • nS'2 • 'i'2'?i • n + 35i • nS'2 • 'i''2'S2 • n + 75i • nS'2 • nvi ■ V2 

- bSi X V2 ■ n §2 X vi ■ ft + Si' (2(S2 x V2y - 3(S2 x n)*t;2 • ft - 4n'S2 x V2 ■ n 



Fig. 3{h) 
Fig. 3(i) 



= -2 



G^m2 



S' 



f (S^' + {S2 X V2y ~ 4n'S°^n^ - 4n'S2 x V2 ■ n 



G'^m2 



Si ■ S2V1 ■ V2 + Si ■ S2V1 ■ nw2 • n — Si • n52 • tTiW2 • n — Si • nS2 • ftvi ■ V2 



+ Si X V2 ■ ft S2 X Hi ■ ft + S 



f {iS2 X rty 



V2 ■ ft + 2n*S2 X V2 ■ ft 



(51) 
nf 

(52) 
(53) 
(54) 



(55) 
(56) 

(57) 



The cubic self-gravitational interaction diagrams all require the application of one-loop tensor integrals, see Ap- 
pendix]^ First, one applies the one-loop tensor integrals, which are required here up to order 4, and then the Fourier 
tensor integrals. Here, we first encounter time derivatives, arising from the spin couplings, propagator correction 
vertices, and time dependent cubic self-gravitational vertices. 



(al) 



6- 



(cl) 



(a2) 



(c2) 



(a3) 



(c3) 



(a4) 



0-x-< 



(a5) 




(d2) 



0- 
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FIG. 4. NNLO spinl-spin2 Feynman diagrams of order G : Cubic self-gravitational interaction. These diagrams should be 
included together with their mirror images. 



The values of the cubic self-gravitational interaction diagrams are then given by 



Fig. 4(al) == 3- 



Si ■ S2 — 25i ■ n S2 ■ n 
S^' ({82 X viY - (82 X V2y + 4(5^2 X nyv2 ■ n - An'S2 x v^ ■ n 



Fig. 4(a2) 



Fig. 4(a3) 



Fig. 4(a4) 



- S^' U(Si X viY - 4(Si X V2y - 16(Si X n)Vi • n + 4(Si x n)V2 • n 
+ 12n*5'i X V2-n) 

35*1 • 52^2 - SS*! • V2S2 ■ V2 - 45*1 • S2{v2 ■ nf + 45*1 • W2S'2 • nv2 ■ n 



SSi X V2 ■ fi S2 X V2 ■ n + S2'' \'i{Si X V2y - 4(5i x n) V2 • n - 8n*5i x V2 ■ n 



&'m2 



S\ ■ S2V1 ■ V2 — Si ■ V2S2 ■ vi — 45*1 X vi ■ n S2 X V2 ■ n 



+ S^' (82 xv2y- An'S2 xv2-n 



G^m2 r 



2r4 



35*1 • S2V2 + 3Si ■ S2V1 ■ V2 + 25*1 • V1S2 ■ V2 — Si ■ V2S2 ■ vi — 3Si ■ V2S2 ■ m 

+ 45*1 • 5*2 {v2 ■ n)^ — 12Si ■ S2V1 ■ nv2 ■ n — %Si ■ V1S2 ■ nv2 ■ n + 451 • V2S2 ■ nvi ■ n 
+ 8Si ■ V2S2 ■ nv2 ■ n — 8Si ■ nS'2 • '!7i'S2 • n — 8Si ■ nS'2 • 'S2'S'i • n + l2Si ■ nS2 ■ W2W2 • n 
— 85*1 • nS'2 • nv\ ■ V2 — 45*1 x V2 ■ n S2 x V2 ■ n — 365'i • nS'2 • n{v2 ■ nf' 



+ 48Si • nS2 • fivi ■ nv2 ■ n 



Fig. 4(a5) 



G^m2 



Si ■ V1S2 ■ V2 + Si ■ V2 S2 ■ vi — 4Si ■ V1S2 ■ nv2 ■ n — 4Si ■ V2S2 ■ nvi ■ n 
ASi ■ nS2 ■ 'S'ii/2 • n — 4Si • nS2 • '?2'?i • n — 4Si • nS2 • nvi ■ V2 



-f 24Si • nS2 • fivi ■ nv2 • n 



Fig. 4(a6) 



2r4 



Si • S2w| - Si • S2W1 • V2 + 6S1 • W1S2 • V2 + 7Si • W2S2 • Vl - 5Si • W2S2 • V2 



(58) 



(59) 



(60) 



(61) 



(62) 



— 4Si • S2 {v2 ■ nf + 4Si • S2W1 • nti2 • n — 24Si ■ V1S2 ■ nv2 ■ ft — 28Si ■ V2S2 ■ nvi ■ n 

+ I6S1 • V2S2 ■ nv2 ■ n — 24Si • nS2 • W1W2 • n — 24Si • nS2 • W2'i''i • n + 12Si • nS2 • W2W2 • n 

4- 8S1 • nS2 • nu| — 24Si • nS2 • nvi ■ V2 — 4Si x W2 • n S2 x V2' n — 36Si • nS2 • n{v2 ■ n) 

+ 144Si • nS2 • nUi ■ nv2 ■ n , 



(63) 
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G'^m2 



Fig. 4(62) = 



2r4 

G^m2 
2r4 



Fig. 4(61) = Z" Si ■ S2V1 ■ V2 ~ Si ■ V2S2 ■ vi — 45i x vi ■ n S2 'x V2 ■ n 

2r** 

+ 5-2* ({Si X viY - in' Si xvi-n)], 

Si ■ S2V1 ■ V2 ~ Si ■ V2S2 ■ vi — 45i X vi ■ n S2 X V2 ■ n 

f(^iS2XV2y~An'S2XV2-n^] , 

Fig. 4(cl) = — S^i • 52^2 - 5Si ■ V2S2 ■ V2 + 8Si ■ S2{v2 ■ nf + 165*1 ■ V2S2 ■ nv2 ■ n 

+ 24:Si ■ nS2 ■ W2t''2 • n + 8Si ■ nS2 ■ nvl — ^Si x V2 ■ n S2 x V2 ■ n 
— 72Si ■ nS2 ■ n{v2 ■ nf' 



(64) 



(65) 



Fig. 4(c2) 



G^m2 



llSi ■ S2V2 — 55*1 • V2S2 ■ V2 — 45*1 • S2{v2 ■ nf' + 451 • ^252 • nv2 ■ n — \2Si ■ nS2 ■ vwl 
165i X V2 ■ ft S2 X V2 ■ n 



Fig. 4(c3) = - 



G^m2 



Si ■ S2V1 ■ V2 + Si ■ V1S2 ■ V2 + 45i • 52^1 • nv2 ■ ft — ASi ■ V2S2 ■ nvi ■ ft 



ASl ■ ftS2 ■ ftvi ■ V2 — 4S'l X Vi ■ ft S2 X V2 ■ ft + ASl X V2 ■ ft S2 x Vi ■ ft 



Fig. 4(dl) = 



G^m2 
2r4 



25i • S2VI - 5i • 52^1 ■ fj2 + Si- V2S2 ■ vi - 25i • V2S2 ■ V2 - 45i • 52(w2 • ftf 



+ 45i • S2V1 ■ fiv2 ■ ft — ASi ■ V2S2 ■ ftvi ■ ft + ASi ■ V2S2 ■ ftv2 ■ ft — 45i x V2 ■ ft S2 x V2 ■ ft 
- 52°' USi X vif - 2{Si X V2f ~ 4(5i x ft^vi ■ ft + 4(5i x n) V2 • ft + ^n'Si x V2 ■ n 



Fig. 4(d2) 



Fig. 4(el) 



G'^m2 
2r4 



Si ■ S2V1 ■ V2 — Si ■ V2S2 ■ vi — 45i xvi-ftS2 X V2 ■ ft 



+ Sf (52 X V2T - An'S2 xv2-n 



G^m2 



Si ■ S2V1 • i/2 - 5i • wi52 • U2 - 25i • U252 • wi + 25i • W252 • ^2 + 45i • 52(1/2 • ftf 

— 45i • 52^1 • fiv2 ■ ft + 45i • wi52 • riU2 • ft + 85i • U252 • ftvi ■ ft — 125i • V2S2 ■ ftv2 ■ ft 
+ 45i • n52 • ■yiW2 ■ ft + 45i • n52 • W2W1 • ft — 85i • n52 • W2W2 ■ ft — 45i • n52 • nw| 
+ 45i • ftS2 ■ ftvi ■ fJ2 — 45i xv2-ftS2 xv2-ft + 245i • n52 • n(w2 • ftf 



245i • n52 ■ ftHi ■ ftv2 ■ ft 



Fig. 4(e2) 
F^g■ 4(/) 



G''m2 



7-4 
G'^m2 



Si ■ S2V1 ■ V2 — Si ■ V2S2 ■ vi — 45i • 52171 • ftv2 ■ ft + 45i • W252 • ftvi ■ ft 



2Si ■ 52^2 + 25i • S2V1 ■ V2 + 35i • V1S2 -^2 + Si- V2S2 ■ vi + 45i • U252 • V2 

85i • S2V1 ■ fiv2 ■ ft ~ 125i • wi52 • ftv2 ■ ft — 45i • V2S2 ■ ftvi ■ ft — 245i • V2S2 ■ ftv2 ■ ft 
125i • n52 • uit/2 • ft — 125i • n52 • fJ2Vi ■ ft — 245i • n52 • 'S'2W2 ■ ft — 125i • n52 • ftv2 
125i • n52 • ftiJi ■ V2 — 85i x V2 ■ ftS2 xv2-ft + 72Si ■ ftS2 ■ ft{v2 ■ ftf 
72Si ■ ftS2 ■ ftvi ■ ftv2 ■ ft 



(66) 



(67) 



(68) 



(69) 



(70) 



(71) 
(72) 



(73) 



C. Order G Feynman diagrams 

For the NNLO spinl-spin2 interaction, we also have 13 diagrams at order G^: three are three-graviton exchange 
diagrams, three are factorizable diagrams with one-loops, and seven are two-loop diagrams, as shown in Fig. [5j All 




.^""\. 



(e) 



0-- 






< 



(fl) 



(f2) 



(gl) 



Q-- 



(di) 



(g2) 



(>-■ 



(d2) 
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(i) 



FIG. 5. NNLO spinl-spin2 Feynman diagrams of order G''\ Diagrams (al), (b), (c), (dl), (d2), (fl), (f2), and (i) should be 
included together with their mirror images. 



of which are new spinl-spin2 diagrams. The three-graviton exchange diagrams, constructed with either one-, two-, 
or three-graviton spin couplings, corresponding to diagrams (al), (a2), and (b), just factorize into a product of three 
Fourier tensor integrals. We then have diagrams (c), (dl), and (d2), which factorize into a product of a one-loop 
diagram and a one-graviton exchange. 

The values of these simpler diagrams are given by 



Fig. 5(al) = 2 



0^17111712 



55i • 5*2 — 135*1 ■ fi S2 ■ n] , 



Fig. 5(a2) = ^^ 
Fig. 5{b) 



I 255i • 5*2 — 575*1 ■ n S2 ■ n] , 
— I 115i • 52 — 275i ■ n b2 ■ n 



13 /^■^ 



( 35*1 • 52 — 55*1 ■ n S2 ■ n 



Fig. 5(c) = _ G\mi + ml) ^ ^^^ ^^ _ ^^^^ . ^ S2 ■ n] , 
Fig. 5(dl) = -S^^^ 



Fig. 5(d2) = -4 



G^mi7n2 



Si ■ S2 — 25i ■ fi S2 ■ n 



(74) 
(75) 
(76) 
(77) 
(78) 
(79) 



As for the two-loop diagrams, with two cubic vertices or one quartic vertex, corresponding to diagrams (e) through 
(i), they contain two- loop Feynman integrals, where we have to deal with three kinds of two- loop integrals: 

1. Factorizable two-loops, as in diagrams (e) and (h), which factorize into a product of two one-loops, such that 
each one-loop can be performed separately. 

2. Nested two-loops, or recursively one-loops, as in diagrams (fl), (12), and (i), on which a one-loop is nested in a 
one-loop, so that they should be performed successively, first the nested one-loop, then the outer one-loop. 

3. Irreducible two-loops, as in diagrams (gl) and (g2), which can be formally reduced, using an integration by parts 
method f50], to a sum of integrals of the two previous kinds, i.e. to a sum of factorizable and nested two- loops. 

The factorizable two-loop diagrams yield here purely short distance contributions, of the form S^'^'{r), which are 
contact interaction terms that cancel out. For all other two-loop diagrams, calculations should be made, keeping the 
dimension d general, and the limit d — > 3 is only taken in the end. For the irreducible two-loop diagrams, which are 
the most complicated to evaluate here, irreducible two-loop tensor integrals of order 4 are encountered, compared 
to the NNLO orbital interaction, i.e. the 2PN potential, which requires only the irreducible two-loop scalar integral. 
These are reduced using the integration by parts method [50] , see Appendix [A] The expressions for the irreducible 
two- loop tensor integrals reductions contain explicit poles in d = 3, see Appendix \^ but these cancel out in the 
dimensional regularization. 



The values of the two- loop diagrams are then given by 
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Fig. 5(e) = 0, 
F^g. 5(/l) = jj g'("^i+^i) ( 3^^ .S2-bS,-nS2-n 



Fig. 5(/2) 



lG3(TOf 



Fig. 5(5l) = 4 
Fig. 5(g2) = 4 



G^TOl 



TO2 



G^mim2 



3Si ■ 5*2 — 55*1 • n ^2 • n ) , 



5i • 52 — 35*1 • n 52 • n) , 



75i • 52 — 205i • n 52 • n J , 



Fig. 5{h) ^ 0, 
i'lf/. 5(ij = — = ( 35i • 52 — 55i ■ n b2 ■ n 



(80) 
(81) 

(82) 

(83) 

(84) 
(85) 

(86) 



IV. SPIN1-SPIN2 ROUTHIAN AND HAMILTONIAN AT FOURTH POST-NEWTONIAN ORDER 

Summing up all of the contributions from the Feynman diagrams, we obtain the NNLO spinl-spin2 interaction 
Routhian. Here, we give the result with the accelerations and precessions dependent terms, and we do not substitute 
in the 5"' dependent terms, using the SSC, so that the 5°* entries are left as independent degrees of freedom, and the 
additional contributions of the field corrections in the 5*^' entries are not taken into account explicitly in this form. 
We also rearrange terms that should be taken with 1 o 2. Our NNLO spinl-spin2 interaction Routhian for a binary 
system of compact spinning objects is then given by 



nNNLO 



G 

8^ 



3Si ■ S2vlvl + 2Si ■ S2(vi ■ V2) 



25i • wi52 • wii^l + 45i • V1S2 ■ W2W1 • V2 - 45i • W252 • wii/i • V2 



— 25i • W252 • V2vi - 95i • S2{vi ■ nyv^ - 95i • 52(1/2 • nYv{ - 245i • 52Wi • 1/21/1 ■ nv2 ■ n 

+ 65i • i7i52 • 1/1(1/2 • n) + 245i • V2S2 ■ wiwi • 1/1/2 • n + 65i • V2S2 ■ 1/2(1/1 • n) + 65i • V1S2 ■ nvi ■ nv2 
+ 65i • V2S2 ■ nv2 ■ nv\ + 65i • nS2 ■ wiwi • iii'2 + 65i • n52 • 1/21/2 • nv\ + 35i • n52 • nvfv2 

— 65i • n52 • n(vi ■ 1/2) — 125i x vi ■ V2S2 x vi ■ nv2 ■ n — 125i x vi ■ V2S2 x 1/2 • nvi ■ n 

— 125i X vi ■ nS2 X V2 • 1/11/2 • n — 125i x 1/2 • n52 x 1/2 • wii/i • n — 125i x vi ■ nS2 x 1/2 • nvi ■ 1/2 
+ 125i X V2 ■ nS2 X vi ■ nvi ■ 1/2 + 455i • 52(1/1 • n)'^{v2 ■ nf' — 305i • V1S2 ■ nvi ■ 7/(1/2 • n)'^ 

— 305i • V2S2 ■ nv2 ■ n{vi ■ nf' — 305i • n52 • wii/i • 1/(1/2 • nf — 305i • n52 • 1/21/2 • 1/(1/1 • nf 

— 155i • n52 • 1/(1/1 ■ n) v^ — 155i • n52 • 1/(1/2 • 11) i": + 605i x v\ ■ r/52 x 1/2 • nv\ ■ r/1/2 • n 

— 605i X 1/2 • r/52 X 1/1 ■ 1/1/1 ■ 1/1/2 ■ ^ + 1055i • nS2 ■ n{vi ■ n) (1/2 • n) 
+ 4 [S^'S^' (4i/i • V2 - 61/1 • r/1/2 • n) 

— S1S2'' ( ^v\v2 — 4i'2i'^ — Qv\n-'v2 ■ fi + Qv\n-' vi ■ n + 611*1;^ 1/2 • n — 6n^V2Vi ■ n + Qn^n-'vi ■ 1/2 ) 

+ (5°* ( (52 X ViYvi ■ V2 - 02 X VlYvl - 02 X l/2)Vi • V2 - v\02 X V2) ■ Vi + vl02 X 1/2) • Vi 



— 3(52 X i/i)Vi • r/1/2 • ^ + 3(52 X i/i)*(i/2 • fi) + 9(52 x i/2)Vi • 1/1/2 • ^ + 3(52 x i/)Vi • r/wj 

— 2)v\S2 X i/i • r/1/2 • 1/ ^ 3i'5^52 X 1/2 • 1/1/2 • 1/ ^ 3i'252 x 1/1 • r/i/i • n + 611252 x vi ■ r/1/2 • 1/ 
+ 3i;252 X V2 ■ nvi ■ n + 3n*52 x 1/2 • i/ii/i • r/ + 3n*52 x 1/2 • i/11/2 • n — 3n*52 x 1/1 • nvi ■ 1/2 
+ 3n*52 X 1/2 • r/1/1 • 1/2 — 15(52 x i/)*i/i • 1/(1/2 • 11) + 15n*52 x vi ■ nvi ■ r/1/2 • n 

— 15n*52 X V2 ■ r/1/1 • 1/1/2 • 11) ) + (1 < — > 2) 
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G_ 

'8^ 



3S'i • 820.1 ■ nv2 + 2Si ■ 820,1 ■ U2W2 • n — 28i ■ V282 ■ V2ai ■ n + 28i ■ 0182 ■ W2W2 ■ n + 28i ■ 1282 ■ ^201 • V2 



— 681 ■ V282 ■ a.iV2 ■ n + 25*1 • V2S2 ■ noi ■ V2 ~ Si ■ 0.182 ■ fiv^ — 5*1 • 715*2 • aiw| — 45i x V2 ■ 0.182 x V2 ■ n 
+ 45i X V2 ■ nS2 X V2 ■ ai ~ 9Si ■ 82 {v2 ■ n)^ai ■ n + 681 ■ V2S2 ■ nv2 ■ nai ■ n + 6Si ■ 1182 ■ W2t^2 • noi ■ ft 
+ 38i ■ 7282 ■ ai{v2 ■ ri-Y + S^i • 0.182 ■ n{v2 ■ nf + 35i • ¥182 ■ nSi ■ nv2 — Q81 ■ 1182 ■ noi ■ W2ti2 • n 

— 158i ■ nS'2 • n{v2 • n) oi ■ n + 128i x n ■ 0182 x V2 ■ nv2 ■ n — 128i x V2 ■ nS'2 x n ■ aiV2 ■ n 

+ 481 ■ 82V2 ■ n — 451 • V282 ■ n — 451 • ^82 ■ V2 + ISS'i • 1182 ■ nv2 ■ n 

+ 4 5°' {{82 X ai)V2 • n - v\S2 x n ■ oi - 11,^82 x 1/2 • Si + 3n^S2 x n ■ aiV2 ■ n) - S"* ({Si x V2yai ■ ft 

— {81 X nySi ■ V2 + 2{Si X diyv2 ■ n + W2>5'i x n ■ oi — n'Si x V2 ■ Si + 2a\8i x V2 ■ n — 3(51 x nyv2 ■ ndi ■ n 
+in'Si X n ■ diV2 ■ nj + 5?' (s^'v2 ■ n + v'^S^^n^ - Sn'S^^v'^ + ^n'S°'n^V2 ■ n ~ 2(^2 x n)' + (^2 x 1/1)^2 • n 

— 3(^2 X V2yv2 ■ n — (^2 X nyv2 + V2S2 X vi ■ ft — V2S2 X V2 ■ n — 11^82 x U2 ■ "Si + 3(^2 x n)*(w2 • n) 



3n*52 X vi ■ nv2 ■ ft + 3n*52 x i!2 ■ fiv2 ■ ft 



2] 



'8^ 



5i • 8201 ■ d.2 — 35i • ai52 • 0.2 + 55i • a252 • Si + 35i • 52ai • fia.2 ■ ft. — 81 ■ 0.182 ■ ^02 ■ ft ^ 81 ■ a.282 ■ ftdi ■ ft 

— Si ■ fi82 ■ did.2 ■ ft— 81 ■ ft82 ■ 0.20.1 ■ ft + 35i • n52 • fldi ■ 0.2 + 45i x ft ■ 0182 x ft ■ 02 — 45i x ft ■ 0.282 x ft ■ di 
+ 35i • ft82 ■ ftdi ■ ftd2 ■ ft 

- 4 [5°* 5°' + Sfn'Sl^n^ 

Ji"* (2(52 X 02)' - (52 X nyd2 ■ ft + n'S2 x n ■ 02) - 8^' (2(5i x ai)* - (5i x Hydi ■ ft + n'Si x rt ■ oA 

-> -> Tfli -> Till -» ^ ^ 

5i X U2 • fli — 5i X n • aiW2 • ft + 2 52 x W2 • Si — 2 — 52 x ft ■ S1W2 ■ ft — 81 x fli ■ {J2 + Si x wi • ftv2 ■ ft 

m-2 2712 



-8 



Gto2 

2r 



-25i xv2-nv2-n + SS'^'v^ + S'^'nVa • rt 



GVI2S1 X S • S2 + [1 



2] 



-\-—Gmim,2 [12S1 • i/iW2 ■ S — 14Si • t/2'i/2 ■ ft + di ■ rtwj — Si • fi{v2 ■ ft) — I2S2 ■ ■i/2'?i ■ ft + I4S2 • fJiVi • n — S2 • ftvi 
8 

+S2 • n(wi • S) ] H — Gmim2r [15Si • S2 — Si • SS2 • S] 



2r4 



115i • 52^2 + 35i • V1S2 ■ fJ2 + 55i • W252 • vi — 165i • W252 • V2 — 225i • 52(w2 • S)^ + 165i • 52i''i • ftv2 ■ ft 



+ 25i • wi52 • Si;2 • ft — 245i • W252 • ftfJi ■ ft + 405i • f]282 • fift2 ■ n + 2Si ■ fi82 ■ vift2 ■ tt — 165i • S52 • f]2'di ■ ft 
+ 165i • n52 • fj2'd2 ■ tt + 7 81 ■ ftS2 ■ ftv2 — 165i • n52 • ftfli ■ V2 — 45i x fJi ■ ft 82 x fl2 ■ ft 
+ lOSi X V2 ■ ft 82 X fji ■ ft — 85i X fj2 ■ ft 82 X fj2 ■ ft ~ 425i • S52 • ft{v2 ■ ft) — 125i • S52 • ftfJi ■ ftv2 ■ ft 
- 45i°*5^' + 16S^'n'S°'n^ + 2Sf (2(52 x viy - 6(52 x ^2)^ - 2(52 x rtyvi ■ n + 11(52 x S)V2 • rt 

-5n'52 X wi • S + 6n*52 x V2 ■ nj - 28^' (2(81 x viY + 2(5i x 1/2)'' - 6(5i x S)Vi • rt - 3(5i x S)V2 • rt 



+5 



— 2n^Si X fji ■ ft + n^ Si x fj2 ■ ft 
G^imf +TO2) 



( 



5i • 52 — 35i ■ ft 82 ■ ft 



^2] 
\ G^mim2 ( 



395i • 52 - 1275i ■ rt 82 ■ rt 



(87) 



Now we can substitute in the lower order EOM, i.e. the accelerations and precessions from the orbital and spin 
interactions, which start contributing at this order, as given in Appendix [BJ This is a common procedure for higher- 
order PN corrections |47h49| . The time derivatives on the temporal spin components, i.e. 5°', are more complicate 
beyond the leading PN contribution from the spin-orbit sector, since they contain contributions from several sectors 
at distinct PN orders, as one would expect from just naively differentiating 5°' ^ Sv. Therefore, in a few terms, 
where they are found, we shifted the time derivatives from them to the other variables, which are physical ones, 
neglecting total time derivatives. Keeping only possible spinl-spin2 terms up to the PN order considered here, and 
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again rearranging terms that should be taken with 1 f^ 2, the expression for the NNLO spinl-spin2 Routhian narrows 
down, and we obtain 



R 



NNLO 
S1S2 



G 

8^ 



35*1 • S2vfv2 + 25*1 • S2ivi ■ V2f - 25*1 • V1S2 ■ viv^ + 45i • V1S2 ■ W2t^i • V2 - 45*1 • V2S2 ■ iTiWi • ^2 



— 25*1 • ^262 • 'V2vf — 9Si ■ S2{vi ■ n)^V2 — 951 • S2{v2 ■ nj'^vf — 245*1 • S2V1 ■ ■!T2Wi • nv2 ■ n 

+ 6Si ■ V1S2 ■ 1^1 (w2 ■ n)'^ + 24S'i • U2'S'2 • wi'S'i • nv2 ■ ft + 6Si ■ V2S2 ■ ■i?2(t?i ■ n)'^ + 66*1 • wiS'2 • nvi ■ nv^ 

+ 6Si ■ V2S2 ■ nv2 ■ nvi + &Si ■ nS'2 • wi^i • nw| + &Si ■ nS2 ■ W2'i'2 • nvi + S^i • nS'2 • nw^w| 

~ 65*1 • nS'2 • n{vi ■ 1/2)^ — 12S'i X vi ■ V2S2 x vi ■ nv2 ■ n — 125*1 x vi ■ V2S2 x V2 ■ nvi ■ n 

~ 125i X vi ■ nS2 X V2 ■ '!7iW2 • n — l2Si x V2 ■ nS'2 x V2 ■ wii/i • n — 12Si x vi ■ nS2 x V2 ■ nvi ■ V2 

+ 125i X V2 ■ nS2 X vi ■ nvi ■ V2 + 45S'i • 52(wi • n)^{v2 ■ nf — 305i • V1S2 ■ nvi ■ n{v2 ■ nf 

— 305i • V2S2 ■ nv2 ■ n{vi ■ nf' — 305i • 7182 ■ 'S'l'iJi • n{v2 ■ nf' — 305i • nS2 ■ '!72'i''2 • ^(^i • n)^ 
~ 155i • nS2 ■ n{vi ■ n)^w| — 155i • n52 • n{v2 ■ n)'^vf + 6OS1 x vi ■ nS'2 x V2 ■ nvi ■ nv2 ■ n 

— 6OS1 X V2 ■ nS2 X vi ■ nvi ■ nv2 ■ n + 105Si • nS2 • n{vi ■ nf'{v2 ■ n)'^ 
+ 4 [Si°*S^'' (4wi • V2 - 6wi • nv2 ■ n) 

— Si Sj'' ( Av\v2 — Av\vi — 6v\n^V2 ■ ft + 6v2n-'vi ■ n + 6n^vlv2 ■ n — Gn'wj^i • n + Qn^n-'vi ■ V2 1 

+ fs?* nS2 X t/i)Vi • V2 - (S2 X Vifvl - (S2 X V2yvi ■ V2 ~ vl{S2 X V2) ■ Vi + V\{S2 X V2) ■ Vi 

— 3(S2 X t/i)Vi • nv2 ■ n + 3(S2 x wi)*(^^2 • nf + 9(S2 x 1/2)^1 • nv2 ■ n + 3(82 x n)Vi • nvi 

— 3U1S2 X vi ■ nv2 ■ n — 3U1S2 X V2 ■ nv2 ■ n — 3t;2S2 x vi ■ nvi ■ n + 6W2S2 x vi ■ nv2 ■ n 
+ 3V2S2 X V2 ■ nvi ■ n + 3n*S2 x V2 ■ wiWi • n + 3n^S2 x -112 • W1W2 • n — 37i*S2 x vi ■ nvi ■ V2 
+ 3n'S2 X V2 ■ fivi ■ V2 — 15(S2 x n)*tTi • n{v2 ■ n) + 15n*S2 x vi ■ nvi ■ nv2 ■ n 

— 157i*S2 X V2 ■ nvi ■nv2-n + 2(S2 x vif - 2(S2 x 1/2)* - 6(S2 x n)''(t/i ■ n - V2 ■ n)\'\ + (1 i — > 2) 

-— -^Si* [i;iW2 • n - Sv^vi ■ n + 2^21/2 • n + n' (vi ■ V2 ~ 



3(wi • ^2 • n - [V2 ■ nf))] + [1 ^^ 2] 



G^m2 
4r^ 



2Si • S2U2 + I6S1 • S2V1 ■ V2 - 6S1 • tTiS2 • V2 + 4Si • V2S2 ■ Vi - llSi • V2S2 ■ V2 



7Si • S2('y2 ■ n) — 42Si • V2S2 ■ nvi ■ ft + 86S1 • U2S2 ■ nv2 ■ n — lOSi • nS2 • ■i/ii72 • n 
I8S1 • nS2 • '!72Wi • n + 38Si • nS2 • W2'i'2 • n + 17 Si ■ 71S2 • nf 2 ^ 30Si • nS2 • nwi • V2 
2OS1 X wi • n S2 X f 2 • n + 6S1 X f2 • n S2 X wi • n + I6S1 x i/ 2 • n S2 x wi • n 
183Si • nS2 • n(w2 • n)^ + 42Si • nS2 • nvi ■ nv2 ■ n 



-2 



9Si°*Sf - 29Si"n*S^^n^' - S?' (4(S2 x viY - 10(S2 x mY - 4(S2 x n)Vi • n + 22(S2 x n)'v2 ■ ■ 
'10n*S2 X fi • n + 5n*S2 x V2 ■ n - {v2 x n)^S2 • n) + S°^ (lO(Si x vi)' - 3(Si x V2)* 
- 21 (Si X n)Vi • n - (Si X ?1)V2 • n - 14n*Si x wi • n + 22n'Si x V2 ■ n - 2(t/i x n)*Si • n 



+5 



+ 4(1/2 X n)'Si • n 
G^imj +m|) 



2] 



Si • S2 — 3Si • n S2 • n 



r 
G^m 



3 STn' + (1 ^^ 2) 

■g — I 55Si • S2 — I9IS1 ■ n S2 ■ n 



2r 



(88) 



From this, one can derive the NNLO spinl-spin2 Haniihonian essentiahy following similar steps to those reported in 
Sec. VI of [25]. These include (1) Legendre transforming with respect to the velocities, (2) substituting in the covariant 
SSC, including the tetrads corresponding to the metric for a binary of spinning black holes in harmonic coordinates, 
up to order G^, and (3) mapping the position and spin from "covariant" to canonical variables, where these mappings 
should be extended to higher orders. In fact, it may be possible to eliminate the S'° degrees of freedom using the 
covariant SSC at the level of the action, and get a reduced Routhian result. In a forthcoming paper, this derivation 
will be shown in full detail, and the resulting Hamiltonian will be compared with the NNLO spinl-spin2 Hamiltonian 
obtained via the ADM canonical formalism in [45| . Alternatively, from this result the EOM can be computed directly, 
using the full spin algebra as in, e.g. |46| , see Appendix [Bl and then like the steps mentioned above for the obtainment 
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of the Hamiltonian one should (1) substitute in the covariant SSC, including the tetrads up to order G^, and (2) map 
the position and spin from covariant to canonical variables to get to canonical EOM. 



V. CONCLUSIONS 

In this paper, we calculated the NNLO spinl-spin2 potential for a binary of compact spinning objects at the 4PN 
order. Such high PN orders are required for the successful detection of gravitational radiation. We have performed 
the calculation using the EFT approach and in terms of the NRG fields. Here, we first demonstrate the ability of the 
EFT approach to go beyond the NLO in PN corrections of spin effects. The NNLO spinl-spin2 interaction sector 
includes contributions from 51 diagrams, of which 45 are pure spinl-spin2 diagrams, while a further 6 arise from other 
sectors, but contribute through the LO spin EOM, that should be taken into account here for the first time. Thus, 
we first encounter here diagrams, that are not pure spinl-spin2 diagrams, but still contribute to the interaction. 39 
new diagrams appeared here, while 6 others contained new ingredients. 

Two features of the spin couplings present the main difficulties. First, the fact that the spin, formally being a tensor, 
is derivative-coupled unlike the scalar mass. This calls for higher-order tensor expressions for all integrals involved in 
the calculations, which add significantly to the complexity of computations. In particular, for the irreducible two-loop 
diagrams, which are the most complicated to compute here, irreducible two-loop tensor integrals up to order 4 are 
encountered. Moreover, the derivative-coupling also allows for time derivatives in the worldline spin couplings, which 
are an additional complication in spin computations. The time derivatives also make the corresponding terms scale 
at higher PN orders. The second obstacle is the fact that the spin couplings contain 5'° entries, which represent the 
redundant unphysical degrees of freedom related with the spin tensor. These are taken as independent degrees of 
freedom throughout the calculation, yet eventually, possibly even after the obtainment of the EOM, they are reduced 
from the final result using some SSC. The 5"* entries also yield contributions of higher PN orders with respect to 
the S^^ spin tensor components. Both features make the PN order of the spin couplings implicit, compared to the 
mass couplings, and the power counting, which is essential in the EFT approach, is more difficult with respect to the 
nonspinning case. 

Moreover, unfortunately, not all of the good attributes of the NRG fields pass on to spin interactions. In particular, 
all possible diagram topologies are realized at each order of G included, as was already illustrated in the NLO spin 
interactions. Also the derivation of the spin couplings is more convenient in terms of the standard Lorentz covariant 
parametrization, and is not so simple and immediate as for the mass couplings. Still, the NRG fields remain advan- 
tageous over other parametrizations, and thus there was no use of automated computations in this work. However, 
for calculations beyond this order with or without spin effects, it is clear that automated computations utilizing the 
NRG fields should be implemented, and would be most powerful and efficient. 

Our final result here can be reduced, and a NNLO spinl-spin2 Hamiltonian may be derived from it. In a forthcoming 
paper, this derivation will be shown in full detail, and the resulting Hamiltonian will be compared with the NNLO 
spinl-spin2 Hamiltonian obtained via the ADM canonical formalism in [i^. Alternatively, the EOM can be computed 
directly from our final result. Future work may shed more light on the equivalence of the ADM canonical formalism 
and the EFT approach in the treatment of spinning bodies, and may result in improvement in both. 
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Appendix A: Dimensional regularization and Feynman integrals 

Throughout the computation of the contributing Feynman diagrams, we encounter in general two types of momen- 
tum integrals that need to be evaluated: Fourier integrals that arise from the Fourier transforms of the propagators, 
and loop integrals, which arise from the cubic and quartic self-gravitational interaction. Both types of integrals 
are evaluated using dimensional regularization [50| . In order to evaluate the Fourier integrals, one should use the 
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d-dimensional master formula for the scalar integral given by 



d^k e''^"- 1 T{d/2-a) fv^Y"^''' 

- ^ ' ' I \ _ (Al) 



(27r)'' (k2)" (47r)''/2 V{a) \ 4 

This formula can easily be derived using Schwinger (or alpha) parameters [50[ . Differentiating the above with respect 
to r, yields the following d-dimensional master fornnilas for the tensor Fourier integrals: 
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Actually, for this paper, the tensor Fourier integrals of orders 5 and 6 were required. We did not include here the 
lengthy expressions for them, containing 26 and 76 generic terms, respectively. 

In order to evaluate the loop integrals, one should use the d-dimensional master formula for one-loop scalar integrals 
given by 
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This formula can easily be derived using Feynman and Schwinger parameters [5Cr]. Similarly, one can also derive the 
following d-dimensional master formulae for the one-loop tensor integrals: 
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In addition, we encounter irreducible two-loop tensor integrals up to order 4. These can be reduced, using an 
integration by parts method [50|, to a sum of factorizable and nested two- loops, as explained in Sec. IIII CI The 
required irreducible two-loop tensor integrals reductions are given by 
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It should be noted that these expressions contain explicit poles in d = 3, but these cancel out in the dimensional 
regularization. 



Appendix B: Accelerations and Precessions 



In the diagrams of Figs. [T] and [2 there appear spin-orbit accelerations Si^go) and precessions S(so)i ^?so)^ ^^'^ 

spinl-spin2 precessions ^(SjSj), 5'?^ g ■,, that should be substituted in the spinl-spin2, and spin-orbit diagram values, 
respectively. LO spinl-spin2 accelerations ^(SiSa) should also be substituted in the orbital interaction diagrams. This 
is in addition to the required substitution of the Newtonian acceleration a/f{\ . We note that the time derivatives on 
the temporal spin components, i.e. S*"*, are more complicate beyond the leading PN contribution from the spin-orbit 
sector, since they contain contributions from several sectors at distinct PN orders, as one would expect from just 
naively differentiating 5°* ~ Sv. Therefore, in a few terms where they are found, we shifted the time derivatives from 
them to the other variables, which are physical ones, neglecting total time derivatives. Thus, the required EOM of 
spin effects are derived here from the LO spin-orbit Routhian in [25|, and the LO spinl-spin2 Routhian in 16]. The 



22 

accelerations are derived from the Euler-Lagrange equations or by applying the variation principle on the Routhian, 
namely, 
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The accelerations required for the NNLO spinl-spin2 potential are then: 
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where we have written them for particle 1, and for particle 2, the equations should just be taken with the particle 

labels exchanged, i.e. 1 o 2. For the last equation here the vanishing of the spin precession S at Newtonian order is 
used. 

The precessions are obtained using Hamilton's equations for the Routhian |46j . namely: 

5 = {5,7^}, S•"* = {5"^7^}, 
where the reduced spin algebra is just obtained from the full spin algebra, e.g. in [46[, such that: 
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The precessions required for the NNLO spinl-spin2 potential are (also given for particle 1, with 1 -(-> 2 for particle 
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where in the last equation for S^^^q-. only the leading PN order contribution is considered here. 
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